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Abstract 



The Kompaneets theory of photon kinetic evolution due to the Compton effect is extended to the case 
of the Vlasov plasma wave oscillations. 

Taking into account Zel'dovich-Levich's approximation we study interaction of accumulating photons 
with plasma in the long wavelength limit. 



Introduction 

To consider the role of the Compton scattering of quanta on non-relativistic electrons, Kompa- 
neets proposed [1] the following kinetic equation for the photon-number distribution function 
n Ko (hco,t) in a unit volume: 



-Jd 3 pJ dW fi Jto (nw,t)(l + n jro (W,t))/ (e) (1) 



dt 

+ fd 3 p I dW n Ko {hu' ',£)(! + n Ko (huj,t))f (e + K{u - to')); 



where n Ko (hu,t) = J d 3 r n Ko (hco, r, t) , and the integrand is the photon density for a given en- 
ergy hu and time t. Here fo{s) is the Maxwell distribution function for temperature T e of free 
electrons, e = p 2 / (2m) ; dW is the differential photon transition probability from one state into 
another due to scattering with electrons. In this equation, emission and absorption processes 
have not been taken into account, so that the transitions are produced exclusively by Compton 
scattering processes. It is easy to check that in stationary case d t n Ko = the solution is the 
Planck distribution: n Ko (M = l/{exp(hu/kT) - 1). 

Kompaneets considered the non relativistic case, i.e., it was assumed that the inequality 
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kT <C mc 2 holds. Introducing a new variable hu/kT = x and notation uo' — uo = A he 
got from (1) in approximation |A| << ui: 

dn Ko ( dn Ko \ h f 3 



+ n Ko (l + n Ko ) j ^Jd 3 P J dWf (s)A 

+ (^ +2(1 + n - ) ^ +n - (1+n - ) ) lwyI d3p I dwfo{£)A2 (2) 



After introduction of dimensionless time parameter t\ = (mc 2 /kT)(l/c)t, where / is the Comp- 
ton range determined by the total cross section 8ire 2 / (3mc 2 ) , and calculation of integrals in (2), 
Kompaneets obtained the following equation: 



dn Ko 1 d 



dt\ x 2 dx 



4 / dn Ko 2 
X \ ~dx Uk ° Uki 
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Later Zel'dovich and Levich [2] found solution of Eq. (3) in the limit n Ko 3> 1 (high-temperature 
regime) and n 2 Ko ^> \dn Ko /dx\, since in this regime Eq. (3) reduces to the inviscid Burgers' 
equation: 

dnzeLe = — — ( x ^ ) (4) 
dh x 2 dx [XUzELj - [) 

They showed that in the absence of absorption the photons undergo a kind of Bose condensation 
in the energy space in the vicinity of zero. This kinetic condensation depends essentially on the 
form of the initial photon distribution. For a certain form of initial distribution, a shock-wave 
as a function of photon energy occurs in the course of its dynamics. The process is extremely 
non-uniform across the frequency spectrum and substantially affected by absorption. Using the 
method of characteristics Zel'dovich and Levich have found solution of Eq. (4) in the following 
form: 

x = F(x 2 n ZeLe ) - 2tix 2 n ZeLe . (5) 

Here F is determined by the initial condition for Eq. (4). According to Eqs. (4)- (5) all points 
on the initial curve F(x) = x 2 n ZeLe (x,t = 0) move along characteristic straight lines parallel 
to the x-axis in the direction of decreasing x with velocity proportional to x 2 n ZeLe . The time 
at which a given point reaches dimensionless energy x = is determined by the expression 
t = F(x 2 n ZeLe ) j \2x 2 n ZeLe ) . Considering a special case of initial condition, which corresponds to 
the Planck distribution n ZeLe (x, 0) = n Ko (x, 0) = l/(exp(:r) — 1) with T ph > T e , one obtains for 
solution (5): 

1 

nzeLe ~ 2^ 

in the frequency region which corresponds to large occupation number, i.e., for x « 1. 




(6) 



Vlasov plasma oscillations 
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After this preliminaries we consider plasma consisting of charged particles with positions r 
and momenta p = raw moving in R 3 . The Vlasov equation [3] for the particle distribution 
function (plasma density) /(r, p,t) is 

d _L + p . d _L _ + g) . = o ( 7 ) 

<9t m dr dr dp 

Here U is an external potential and U is an electrostatic self-consistent plasma potential which 
is defined by U(t,r) = (N /V) J $(r — r')f(r',p',t)d 3 r'd 3 p', where $(r) = e/r. Let us present 
the eventual solution of (7) in the following perturbational form 

uf* 

/(r, p,t) = / (p) + ///* (r, p,t), ^- « 1, (8) 

Jo 

where /o(p) is the space homogeneous Maxwell distribution for electrons. Then the fluctuating 
electrostatic field can be written as 

E M) = ~J ^^r(r',p',t)rfVrfV. (9) 

From Eq. (9) we see that E and /* are of the same order. In what follows, we use the velocities 
v = p/m instead of the momenta. Then using the form of Eq. (8) we get the linearized Vlasov 
equation of the order 0(/x) 

df* df* eE df 

+ v 'l^ 77~ = ' 10 

at or m ov 

divE = -Anej f*(r,v,t)d 3 v. (11) 

Here 

To be specific consider plasma oscillations propagating in r 1 -direction: 

f*(r,v,t) = fkuWexp^iwt + ikri), (12) 
E Tl (r,t) = E kuJ exp(-iut + ikr^, 
E X2 = E r3 = 0. 

In particular we are interested the case of the long wave-length limit for A = 2ir/k. Then 
following Kvasnikov [4], we get from Eqs. (12) in this limit (k — > 0): 

g 

/*( v ^) = /ow(v) exp(-ujt) = y(*)ui/o(v), <p(t) = --^f smuj t. (13) 

Here u is Langmuir frequency. For the photon density we look for the representation which 
is a perturbation of the Kompaneets solution by the Vlasov oscillations: 

n(x,ti) — n Ko (x,ti) + fj,n*(x,ti) with n*(x,0) — 0. (14) 
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Here n Ko (x,ti) is the solution of Kompaneets equation (3). To proceed further we return to 
Eq. (1) for the distribution function of photons, but instead of the Maxwell function distribution 
for the electrons, we consider the perturbed solution (8), where /* is the solution of Vlasov 
equation (10). Then we obtain for unknown function n*(x, t\) , the following linearized equation: 

dn* d 2 n* . n . 1 . . dn* 

Wi = + (7i + 272^ + 1)) 

+ (^ 1 (l + 2n Ko )+ l2 ^l + 2n Ko + 2^) ) n* +<£>(v 1 ,x,t 1 ), (15) 

where 




$(«i,a;,fi) = ^ + n Jfo (l + n Jfo )J — J 'd 3 p J 'dW<p(t)v 1 f (v)A + 



and 



Ti = I d3 " I dW M v ~l A = x< - 4 ~ x )> 



Equation (15) is a linear parabolic equation with non-constant coefficients. Nevertheless, one 
can find solutions in some special cases. Consider them only for some partial case when the 
"shift" term dn*/dx and the "force" which proportional to n* are much smaller than the 
diffusion term. Then 



dn* 9 d 2 n* 



dtl = x ^ + Hvi ' xM) - (17) 

In the framework of the Zel'dovich - Levich approximation n Ko ^> 1 and n\ o ^> \dn Ko /dx\, it 
is straightforward to find that 

<S>(vi,x,ti) = —Lp(ti)xn 2 ZeLe (x,ti), where n ZeLe (x,ti) has the form (6). To find solution of 
Eq. (17), let us change variables: z = —lax, r — t\, and n*(z,r) = n*{x,t\). Then the 
equation for n*(z,r) takes the form 

dn* d 2 n* dn* 

— — = — — + — \-$(vi,z,t), -oo < z < +oo, < T < +oo 

OT OZ A OZ 

n*(z,0)=0 (18) 

After the second change of variables: n*(z,r) = exp(— z/2 — t/A)v(z, t) , we get for v{z,t) 
diffusion equation: 

dv d 2 v 

^ = ^2+ ^ Z > T ) eX P(^/ 2 + T / 4 )' "( Z > ) = - ( 19 ) 



4 



The solution of Eq. (19) is 

A 



n [z,t) 



(----) 

2 4>- 



rT poo 

— I dr' dz' sin u T'n 2 ZeLe (z' , r') x 

TC JO J-oo 



2v^F 

e (¥+|) e 3^y( v / r _ r /)-i (20) 

where A = el / {u^kTc) . 

To calculate integral (20) we consider n ZeLe (x,ti) in the following two limit cases (0 < x < 1): 

• The first one is the case of small times t\ < 1/2 + x/4. Then it is straightforward to 
show that 4t 1 x/(2t 1 - 1 - §) 2 < 1, and TizeLe{x 7 ti) — x(i+x/2—2ti) ' For simplicity let 
us consider the Zel'dovoch-Levich approximation, when n ZeLe (x,ti) = ^jz^Tj ■ Then one 
obtains $(^i,x,ti) = A sinc^i x(l ^ 2ti)2 ■ 

Remark: In this approximation the last formula for n ZeLe manifest a singularity at 
t\ = 1/2. Zel'dovich and Levich have mentioned that a kind of photon Bose conden- 
sation is occurring just at this moment, see also [5], [6]. They have found that it is a 
minimal critical time in which the photon state with x = starts filling up. 
For time t± < 1/2 one obtains from (20) the solution of Eq. (17): 

A „ , 0i N f 1 sin(f - <frj) 
-2ti rj 



n*(x,t 1 ) = ^-(l + 2t 1 )/ 1 Sinr2 ° ^ V) d V 
2x Ji-2ti n 



. UJ [l 1 COS U} T]' , / 

sm — / ■ — dri 

2 i^(i-2ti) rf 

ojq sincuoV 



f— sinw ?)' , 
cos — / — dri 

2 iif(l-2*i) 77' 

For the second limit case, t >> 1 , one gets n ZeLc (x, ti) = 2/a; 2 with $(i>i, x, ii) = 4A smcu Q ti . 
Then 

^_ e/ (12sino;oti - a; cosa; ti) _ ftu; _ mc 2 Z 
" lX,tlj "a;oA:T ^(144 + o; 2 ) > x — kT kT c 



We notice that, in contrast to the first case, the behavior the photon distribution when x goes 
to zero is different from the Planck distribution. It is in fact a dominant contribution as x 
approaches zero. 



Conclusion 

We found the photon function distribution perturbed by the Compton scattering in plasma is: 

n(x, = n ZeLe (x, ti) + fin*(x, ti). (22) 

We have shown that the solution exhibited by Zel'dovitch-Levich might be unstable at suffi- 
ciently large times around x = where photon Bose condensation occurs. This instability is 
generated by the diffusion term in the equation for the perturbation of the photon density by 
plasma wave oscillations. 
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